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ABSTRACT 

We study an assembly-type bias parametrized by the dimcnsionless spin parameter 
that affects massive structures. In numerical simulations higher spin haloes are more 
strongly clustered than lower spin haloes of equal mass. We detect a difference of over 
a 30 per cent in the clustering strength for dark matter haloes of 10 13 - 10 14 M Q , 
fH | which is similar to the result of Bett et al. We explore whether the dependence of 

clustering strength on halo spin is removed if we apply the redefinition of overdensity 
peak height proposed by Lacerna & Padilla (Paper I) obtained using assembly ages. 
We find that this is not the case due to two reasons. Firstly, only a few objects of 
low-virial mass are moved into the mass range where the spin introduces an assembly 
bias after using this redefinition. Secondly, this formalism does not alter the mass of 
massive objects. In other words, the sample of haloes with redefined mass M in the 
high-mass regime is practically the same as before the redefinition of peak height, and 
thus the clustering behaviour is the same. We then repeat the process of finding the 
redefined peak height of Paper I but using the spin. In this case, the new masses show 
no spin-related assembly bias but they introduce a previously absent assembly bias 
with respect to relative age. From this result, we conclude that the assembly- type bias 
with respect to the halo spin has a different origin than with respect to assembly age. 
The former may be due to the material from filaments, which is accreted by massive 
haloes, that is enhanced in high-density environments, thus causing more extreme spin 
values without significantly changing the formation age of the halo. In addition, the 
estimates of the mass of collapsed structures in numerical simulations could be lower 
than the true mass, even in cluster-size haloes. High-mass objects may correspond, 
in some cases, to a different peak height than that suggested by their virial mass, 
providing a possible explanation for the assembly bias with respect to spin. 
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1 INTRODUCTION 

The new generation of numerical simulations of high resolu- 
tion have shown that the large-scale clustering of haloes of 
a given mass varies significantly with their assembly history 
(Gao & White 2007). This effect, which is not expected from 
the extended Press— Schechter theory (Bond et al. 1991), 
was termed 'assembly bias', and is found when measuring 
the amplitude of clustering as a function of halo properties 
such as formation time and halo spin at fixed halo mass. 

In the first paper of this series, we presented a new ap- 
proach to estimate the overdensity peak height with the aim 
to understand the assembly bias effect (Lacerna & Padilla 
2011, hereafter Paper I). The method consisted in redefin- 
ing the overdensity that characterizes each object using the 
information of its virial mass and relative age. This new 



definition is proposed as a better alternative than the virial 
mass as a proxy for equivalent peak height. 

In this letter we will investigate the prevalence of the as- 
sembly bias with the proposed redefinition of an overdensity 
peak height of Paper I using the dimensionless spin param- 
eter. Works by Bett et al. (2007, hereafter B07) and Fal- 
tenbacher & White (2010) have detected an assembly-type 
bias in massive dark matter (sub)haloes using this property. 
They found that higher spin objects are more strongly clus- 
tered than lower spin objects of equal mass. 

The outline of this work is as follows. Section[2]describes 
the data, and the dimensionless spin parameter is defined in 
Section [3] The detection of the assembly bias in our simu- 
lation using halo spin is shown in Section 3] and Section [S] 
shows the results using the overdensity peak height proxy. 
The discussion and conclusions are presented in Section [6] 
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2 DATA 

The numerical simulation we use in this letter is called 
STAND, and consists of a periodic box of 150 h~ x Mpc on 
a side. It contains 640 3 dark matter particles with a mass 
resolution of ~10 9 h -1 Mq. The cosmology used in this sim- 
ulation is Q. m — 0.28 (with a baryon fraction of 0.164), Qa 
— 0.72, as — 0.81, h — 0.7, and spectral index n s = 0.96. 
The chosen value of as is in good agreement with recent 
results of WMAP (Komatsu et al. 2011). The simulation is 
run from redshift z = 73.5 using the initial conditions of 
GRAFIC2 (Bertschinger 2001) and the public version of the 
GADGET-2 code (Springel 2005). 

In 99 snapshots of the simulation, dark matter haloes 
are identified as structures that contain at least 10 parti- 
cles using a friends-of-friends (FOF) algorithm (Davis et al. 
1985). Another algorithm (SUBFIND; Springel et al. 2001) 
is applied to these groups in order to find substructures with 
at least 10 particles. 

The catalogues of haloes and subhaloes are used to con- 
struct merger histories, over which the SAG2 model by La- 
gos, Cora, & Padilla (2008; see also Lagos, Padilla, & Cora 
2009) is run to produce a galaxy population. In this letter, 
we only use the central galaxy of the most massive subhalo 
within a FOF halo. 



3 DIMENSIONLESS SPIN PARAMETER 

The dimensionless spin parameter A is typically defined as 
(Peebles 1971) 



A = 



J\E\ 



1/2 



-1 



< 0.4 
I 0.2 
M 





i 


— o-- 


| ' 

- 50% 


highest spin 






B-- 


50% 


lowest spin 

- 


_ 






'x's. 


- 




STAND 






H) 


10 13 


- 10 14 








: — 1 1 h 


— 1 1— 

B07 

- i-p _ 


H 1 1 1 


h 

STAND 


— 1 1 1 — : 
























1 . 





0.5 



1 



1.5 



log r/h _I Mpc 



GM, 



5/2 



(1) 



Figure 1. Main panel (top): two-point cross-correlation function 
for haloes in the STAND simulation. The result for the 50 per 
cent highest spin haloes is represented as dot-dashed red lines, 
whereas that for the 50 per cent lowest spin haloes appears as 
dotted blue lines. Error bars were calculated using the jackknife 
method. Lower panel (bottom): ratio between the bias of high and 
low spin objects in the STAND simulation (solid lines) and in B07 
(dashed lines). At large scales, r > 3 h~ 1 Mpc, both simulations 
show a higher clustering for high spin haloes with respect to low 
spin ones with a similar signal, which is over 30 per cent. 



where J is the magnitude of the angular momentum vector, 
E is the total energy, G is the gravitational constant, and 
Mh is the mass of the FOF halo defined by the number of 
dark matter particles, N p . The angular momentum is 



T M h ^ 



(2) 



where r; and Vj are the position and velocity vectors of 
particle i relative to the halo centre of mass, respectively. 
Furthermore, the total energy of the system is given by E — 
T + U, where T corresponds to the kinetic energy and U to 
the potential energy. 

We select haloes with number of particles N p ^ 300 be- 
cause they have reliable estimates of this parameter (B07). 
The number of selected objects using this condition is 31,449 
haloes, which corresponds to four per cent of the total sam- 
ple. As noted by B07, there is an important scatter at higher 
spins even for N p 300 due to small objects with high ve- 
locity dispersions, caused by the proximity of more-massive 
haloes. This velocity contamination increases the kinetic en- 
ergy of small haloes, T, resulting in larger values of spin. To 
overcome this problem, they suggested the quasi-equilibrium 
criterion to remove anomalous spin haloes. This is a cut of 
the instantaneous virial ratio of halo energies 2T/U + 1. 
Typically, a virialized object has a value around zero, and a 
gravitationally bound object has value ^ — 1. We select the 
haloes according to 



^ 2T 



(3) 



where the limit value is Q — 0.5 (B07). This allows us to 
remove objects with anomalous spins. The final sample con- 
tains 29,633 haloes. 

Although the spin A has little dependence on mass, if 
any, we define the relative spin parameter, in analogy to the 
relative age parameter 5t of Paper I, as 



A, - (A(A/)) 



(4) 



where, for the i th object, Xi is its dimensionless spin, 
(A(M)) is the median spin as a function of FOF host halo 
mass M, and a\(M) is the dispersion around the median. 
Then, positive values of 5\ correspond to high spin objects, 
whereas negative values of 5\ are related to low spin objects. 



4 ASSEMBLY BIAS USING HALO SPIN 

In this section, we will use the two-point correlation function 
to test whether dark matter (DM) haloes from the STAND 
simulation show an assembly-type bias with the halo spin. 
Figure [1] shows the two-point cross-correlation function 
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of DM haloes of equal mass but different relative spin pa- 
rameter given by equation Q. The centres to perform £(r) 
are haloes selected in the mass range 10 13 - 10 14 hr 1 M , 
and the tracers are all the haloes in the final sample (see 
Section [3} in order to obtain a better signal-to-noise. The 
50 per cent high spin haloes are more strongly clustered at 
large scales than the 50 per cent low spin haloes, thus show- 
ing an assembly-type bias. The lower panel shows the ratio 
between the bias of high and low spin objects in the STAND 
simulation (solid line). This ratio is calculated as 



JH .high x 
bH,low x 



H H' ,highx 



HH'low. 



(5) 



where the subscript H refers to centres and H' to tracers. 
We compare to the autocorrelation function results of B07 
shown as a dashed line in the lower panel. In this case, the 
ratio between the bias of 50 per cent high and 50 per cent 
low spin objects is calculated as 



bn,h 



OH,low x 



£,HH,hi(lh x 



(6) 



As can be seen from the lower panel, at scales r > 3 
h~ Mpc, both simulations show a higher clustering for high 
spin haloes with respect to the low spin ones (over a 30 per 
cent difference) with a comparable signal. 



5 OVERDENSITY PEAK HEIGHT PROXY 

In Paper I, we presented an approach to trace the assembly 
bias effect. This consisted in measuring the total mass in- 
side spheres of different radii around semi-analytic galaxies, 
which in some cases returned larger values than the virial 
mass of the host dark matter halo, using two free parameters 
to introduce a dependence of this radius on mass and age. 
This model defined a new overdensity peak height for which 
the large-scale clustering of objects of a given mass did not 
depend on the age. The parametrization helps to character- 
ize the environment and to understand what is behind the 
assembly bias (see Paper I for more details). 



5.1 Parametrization using the relative age 

We will first use the same approach of Paper I in order to 
find the best-fitting parameters so that £ (r) does not show a 
dependence of the large-scale clustering on age. We will then 
use the parametrization to study whether the dependence on 
spin is still present. 

For this purpose, we match up the position of cen- 
tral semi-analytic galaxies with the position of the selected 
haloes in order to estimate the relative age of these haloes 
using 



U - (t(M)) 
at{M) ' 



(7) 



where, for the i th galaxy, ti is its mass-weighted stellar age, 
(t(M)) is the median stellar age as a function of host halo 
mass M, and a t (M) the dispersion around the median in 
units of time. The procedure is successful in the 99 per cent 
of the cases, so that the final sample contains 29,376 haloes 
for which the results shown in Figure [1] are the same. In 
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Figure 2. Two-point cross-correlation function for haloes of equal 
redefined mass but different spin and different age from the 
STAND simulation. The best-fitting parameters are a = —0.36 
and b = —1.52, where the reduced x 2 statistics is performed us- 
ing the spin (Section 15.21 1. Error bars are calculated using the 
jackknife method. The high-spin population is represented as a 
dot-dashed green line, whereas the low-spin population is repre- 
sented as a dotted cyan line. The old population is represented 
as a short-dashed red line, whereas the young one appears as a 
long-dashed blue line. For clarity, these two lines are moved -0.5 
dex in the vertical direction. Lower panel: ratio between the bias 
of high-spin and low-spin objects (solid line) and between the bias 
of old and young objects (dashed line). The dotted line represents 
the ratio equal to unity. 



the other 257 objects, which include some massive haloes 
of 10 14 ' 5 fc _1 tfg, the position of the central galaxy is very 
far away from the central postion of the halo. We use the 
condition d ^ 250 h~ x kpc, where d is the distance between 
the position of the central galaxy and centre of mass of the 
DM halo. We do not assign the age of the central galaxy 
to its host halo at distances beyond this limit. Again, these 
'failures' correspond to just one per cent of the haloes. The 
resulting stellar age as a function of halo mass is qualita- 
tively similar to that shown in Paper I. On average, more 
massive objects have older stellar ages. 

We redefine the overdensity peak height within spheres 
around haloes of mass Mh, for which the radius of a sphere 
(in h~ Mpc units) is parametrized as 



r = a S t + b log 



M h 
M nl 



(8) 



where M n i is the non- linear mass defined by Seljak & War- 
ren (2004), log(M ni //i _1 M ) = 13.24 for the choice of cos- 
mological parameters in the STAND simulation. The free 
parameters are a and b. The new peak height proxy will be 
the mass M enclosed within this radius. It is assumed that 
if r is smaller than the virial radius or if M is smaller than 
the halo mass, then M = Mh- The \ 2 statistics to find the 
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best-fitting set of free parameters is denned as 

N 



Xf(r) 



N 



(9) 



where £hi g h(r) and ^ offi (r) is the result for old and young 
haloes, respectively. The error is a 2 



'thigh +CT «(o 

the first term is the jackknife error for £hi g h(r) and the 
second for fi ow (r). The symbol ridof denotes the number 
of degrees of freedom (ridof = 6). This statistics is calcu- 
lated within 1 ^ r/h~ x Mpc ^ 13 over three mass bins (the 
first, second and third terciles of the mass distribution), i.e. 
N = 3. 

The best-fitting values are a — 0.08 and b = —0.37. 
As in Paper I, the assembly bias with respect to age is not 
present at large scales after performing this procedure. How- 
ever, the differences in the clustering strength of haloes of 
equal high mass but different spin will still be present be- 
cause just one halo of low virial mass is considered in this 
high-mass range after using the parametrization. Also, as 
mentioned in Paper I, this formalism does not alter the mass 
of objects with high mass. Therefore, the sample of haloes 
with masses between 10 13 and 10 14 h _1 MQ is practically the 
same, thus showing the same clustering behaviour as before. 

5.2 Parametrization using spin 

To trace the assembly bias with the dimensionless spin pa- 
rameter, we use the reduced \ 2 statistics defined in eq. ((9}, 
but in this case £hi g h(r) and ft OI1I (r) correspond to the cross- 
correlation function for high-spin and low-spin haloes, re- 
spectively. This statistics is calculated within 1 sC r/h" 1 
Mpc ^ 13 over two mass bins, 10 12 - 10 13 h~ 1 M Q and 
10 13 - 10 14 hT x M B , i.e. N = 2. 

The best-fitting values are a = —0.36 and b = —1.52. 
After applying the redefinition of mass there is no difference 
in the clustering amplitude for massive haloes of high and 
low spin at large scales, as can be seen in Figure O In con- 
trast, it can be seen in the same figure an artificial assembly 
bias with respect to age using this set of parameters. Note 
that the assembly bias for objects of equal mass but different 
age affects low- mass haloes (Gao, Springel & White 2005). 
Therefore, the formalism applied to the assembly-type bias 
with the dimensionless spin parameter A is not related with 
the successful overdensity peak height defined in Paper I. For 
example, Figure [3] shows that the mean radius for haloes of 
intermediate mass is 10r„i r (solid line), which is larger com- 
pared to 1 - 4r V i r described in Paper I. Furthermore, a few 
very massive objects changed their masses (dashed line), an 
aspect that was not found in the redefinition of the overden- 
sity peak height that explains the assembly bias using the 
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Figure 3. Distribution of radius r of eq. J8} in units of viral radius 
r v i r for two ranges of halo mass Mf, from the STAND simulation 
for the parametrization using spin (see Section 15. 21 1. Haloes of 
intermediate mass (solid line) have mean radius of 10r„; r , which 
is larger compared to 1 - 4r„; r described in Paper I. Also, a few 
very massive objects (dashed line) have r > r v i r , in contrast to 
the redefinition of the overdensity peak height that explains the 
assembly bias using the age. 
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6 DISCUSSION AND CONCLUSIONS 

Perhaps the assembly-type bias found in massive objects of 
equal mass but different spin values is related with the mas- 
sive neighbours responsible for the truncation of the growth 
of smaller objects reported in Paper I. High mass haloes 
reach virial masses according to theoretical expectations, 
but those in high density environment accrete an impor- 
tant fraction of material from filaments (Hahn et al. 2009). 



Figure 4. Median overdensity as a function of relative spin 
parameter, 5^, for massive haloes in the STAND simulation. 
Open circles (dashed line) correspond to normal &\ values. Filled 
squares (solid line) correspond to randomly swapped 8\ values 
(see details in Section [6}. The error bars correspond to the error 
of the median in both cases. 
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This highly collimated material could cause changes in the 
angular momentum of these haloes affecting their spin. On 
the other hand, low-mass objects of equal mass do not show 
important differences in clustering when they are split in 
samples of high and low spin A showing a lack of dependence 
with the environment (B07). However, we have already men- 
tioned that small haloes in the vicinities of massive neigh- 
bours are discarded in this kind of analysis since they suffer 
velocity contamination that increases their kinetic energy, 
thus increasing their spin values. Probably, massive haloes 
fed by material in the region of a truncated smaller halo 
suffer a similar contamination that alters their spins. 

The environmental dependence of massive (f3 ^ 
log(M h //i -1 M Q ) < 14) haloes is shown in Figure g| This 
dependence is estimated as the median overdensity (p ol , er ) 
as a function of relative spin parameter, 8\. The density is 
calculated using all the available haloes in the simulation. 
We then estimate the overdensity with respect to the me- 
dian. As can be seen, a smooth dependence with environ- 
ment is found for low-spin and high-spin objects of equal 
mass (dashed line). It is worth to mention that extreme val- 
ues of spin are represented by higher or lower values of S\. 
Therefore, it is plausible that the assembly-type bias found 
by using the spin parameter is due to a density enhancement 
that affects the spin values of massive objects. 

In order to remove this dependence with environment, 
the relative spin parameter 8\ is randomly swapped among 
haloes in mass ranges of 0.1 dex (see the resulting solid line 
in Figure [4]), without altering their positions. After perform- 
ing this procedure, the halo clustering dependence on spin is 
practically absent, as can be seen in Figure [5] The average 
difference in clustering strength between the fifty per cent 
highest spin haloes and the fifty per cent lowest spin haloes 
is smaller than a 10 per cent, showing that the density en- 
hancement can affect the properties of massive objects such 
as the spin. 

Therefore, in the context of the peak formalism, the 
large amount of material from filaments from the cosmic web 
that a massive halo accretes in high density regions might 
also affect its final mass, thus introducing a spurious appar- 
ent change in the bias. Furthermore, tools used to identify 
objects in numerical simulations, such as FOF or SUBFIND, 
often define structures located very close to a massive object 
as a separated identity. Anderhalden & Diemand (2011) de- 
veloped a formalism to identify lost particles, where many 
small structures actually belong to a neighbour massive ob- 
ject. As a consequence, the mass estimation of structures in 
numerical simulations, even in cluster-size haloes, could be 
lower than the actual mass. It is remarkable that the age of 
massive objects is not strongly affected by this discrepancy 
in mass; however this phenomenon is important in proper- 
ties that are directly related with groups of DM particles 
in high density environments such as the spin, which can 
produce an assembly bias effect. 
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Figure 5. £(r) for massive haloes from the STAND simulation. 
Relative spin parameter 8\ is randomly swapped among haloes 
in mass ranges of 0.1 dex. The result for 50 per cent highest spin 
haloes is represented as dot-dashed red lines, whereas that for 50 
per cent lowest spin haloes appears as dotted blue lines. Error 
bars are calculated using the jackknife method. Lower box: ra- 
tio between the bias of high and low spin objects in the STAND 
simulation (solid line) and in B07 (dashed line). The average dif- 
ference in clustering amplitude is smaller than a 10 per cent when 
halo spins of the STAND simulation are randomly distributed. 



for this work were performed using the Geryon cluster at 
the Centre for Astro- Engineering at UC. 
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